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Non-linear Zeeman splitting of neutral excitons is observed in composition engineered In^Gai— x As 
self-assembled quantum dots and its microscopic origin is explained. Eight-band k • p simula- 
tions, performed using realistic dot parameters extracted from cross-sectional scanning tunneling 
microscopy, reveal that a quadratic contribution to the Zeeman energy originates from a spin depen- 
dent mixing of heavy and light hole orbital states in the dot. The dilute In-composition (x < 0.35) 
and large lateral size (40 — 50 nm) of the quantum dots investigated is shown to strongly enhance 
the non- linear excitonic Zeeman gap, providing a blueprint to enhance such magnetic non-linearities 
via growth engineering. 



Over the past decade, semiconductor quantum dots 
(QDs) have attracted significant interest, mainly due to 
the prospect for their use as integrated, electro-optically 
addressable quantum systems capable of storing and 
processing quantum information. 1 Quantum information 
processing requires the possibility for selective manipu- 
lation of a specific spin qubit within a quantum register. 
Such selective addressing using conventionally applied 
techniques, such as electron spin resonance,^ is rather 
challenging and can more conveniently be achieved with 
recently proposed electrical methods for spin control via 
Lande g-tensor modulationP"^ These approaches exploit 
tuning of the magnetic response by pushing the carrier 
envelope function into different regions of composition- 
engineered quantum dot nanostructures. Although elec- 
trical g-factor modulation has been successfully demon- 
strated in composition-graded AlGaAs quantum wells 7 
and vertically coupled InGaAs QD- molecules^ for sev- 
eral years, significant tuning of g-factors in indi vidual 
self-assembled QDs was achieved only recentlyP^^ It was 
shown that a static electric field can be applied to quench 
the orbital angular momentum of the confined carriers, 
thus, modifying the g-factor!^^ Since most of the ex- 
periments addressing the spin of confined carriers are per- 
formed in magnetic fields it is important to develop a mi- 
croscopic understanding of how magnetic fields influence 
the quantum properties of the orbital states. 

In this paper, we report strong magnetic field induced 
tuning of the exciton g-factor in composition-engineered 
In^Gai-^As-GaAs self-assembled QDs. By comparing 
our experimental results with realistic eight-band k • p 
simulations performed using QD size, shape and compo- 
sitional information obtained from cross-sectional scan- 
ning tunneling microscopy (X-STM) we identify the ori- 
gins of the magnetic field-dependence of the g-factor. Our 
results show that magnetic fields influence the excitonic 
g-factor via a mechanism that differs fundamentally from 
the case of static electric fieldsP^In particular, the com- 
bination of the dilute In-composition (x < 0.35) and 



comparatively large lateral size of the QDs leads to spin 
selective mixing of the lowest energy heavy hole (HH) 
and light hole (LH) orbital states, the strength of which 
is controlled by the external magnetic field. This gives 
rise to a quadratic Zeeman spin-splitting - a phenomenon 
previously observ ed on ly in semiconductor quantum wells 
and sup er lat t ices ! 1 3 1 1 4 1 
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FIG. 1. (color online), (left panel) Topography X-STM 
image of a conventionally grown quantum dot and a dot 
grown using identical conditions but including the In-flush 
technique, (right panel) Outward relaxation analysis of the 
two quantum dots from which the absolute In-concentration 
and spatial composition profile in the dots was extracted. 

The samples investigated were electrically tunable 
GaAs n-i-Schottky photodiode structures into which 
two different types of QDs were embedded in the i- 
region. Both samples were produced using molecular 
beam epitaxy under nominally identical growth condi- 
tions. The first sample was grown utilizing the con- 
ventional Stranski-Krastanov growth, whilst in the sec- 
ond was utilized the partially covered island (PCI) "In- 
flushing" method. 15 In the following, the sample con- 
taining the conventionally grown dots is referred to as 
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non-flushed and, analogously, the sample containing the 
dots grown with the PCI technique as flushed. Both 
samples had a single layer of In^Gai-^As self-assembled 
QDs grown in the 140 nm thick z-region, with a relatively 
high growth temperature of 590 °C. The QD layer con- 
sisted of 8 ML of In^Gai-^As with a nominal In-content 
of x = 0.50, deposited at a rate of 0.41 ML/s and an 
As overpressure of 1.5 x 10 -5 mbar. The comparatively 
high growth temperature is expected to lead to an aver- 
age In-content that is lower than the nominal x = 0.50, 
due to the combined effects of In-desorption, 16 interdiffu- 
sion with the GaAs matrix material and In segregation.^ 
Comparison of our results with the simulations provides 
strong support for this expectation, showing that the 
strong tunability of g ex is inextricably linked to a low 
average x ~ 0.35 and the In-Ga alloy profile. For the 
flushed sample a growth interruption was included after 
the QDs had been partially capped with a 6 nm thick 
GaAs layer. During this growth interruption the tem- 
perature was increased to 650 °C and kept constant for 
30 s. After this, the temperature was again lowered to the 
nominal growth temperature and an additional capping 
layer of GaAs was deposited. X-STM measurements 18 
revealed that the QDs from the high density regions of 
the wafers exhibit inhomogeneous In-composition profiles 
with a relatively large cross sectional size of 40 — 50 nm 
and a height of 4 — 8 nm. Typical topography X-STM im- 
age of two representative QDs is shown in the left panel 
of Fig. [I] Two types of QD composition profiles were 
successfully fitted to the measured outward relaxation of 
the cleaved facet: (i) linearly increasing In-concentration 
from the base to the apex of the dot and (ii) an inverted 
trumpet like In-distribution.^HIII The results from the 
later are shown in the right panel of Fig. [I] Here it 
should be noted that X-STM outward relaxation analy- 
sis can yield approximately similar concentration profiles 
that match to the surface relaxation and that the method 
itself does not provide an answer to which exact com- 
position profile applies to the QDs studied. 22 However, 
our experimental observations were found to be in good 
accord with theory only using the inverted trumpet In- 
distribution profile. The outward relaxation simulations 
revealed that for the conventionally grown non-flushed 
QD shown in the left panel of Fig. [I] the In-concentration 
in the apex is x apex = 0.35, reducing to x m%n = 0.22 
at the base. In contrast, for the flushed sample the In- 
concentration at the apex is slightly lower [x avex = 0.30) 
due to desorption of In during the flush step. 

Optical characterization of the quantum dots was per- 
formed at low temperatures (10 K) using a confocal mi- 
croscope. The microscope was placed in a superconduct- 
ing magnet enabling application of magnetic fields up to 
15 T in Faraday configuration. Typical photolumines- 
cence (PL) spectra are presented in Fig. [5] from the neu- 
tral exciton X° of three representative QDs from the non- 
flushed sample. Many quantum dots (> 30) were studied 
and the results summarized in Fig. [2] illustrate the full 
range of behaviors observed. The polarization-resolved 



PL spectra from the different QDs labeled as QDa, QDb 
and QDc reveal a substantially different behavior of the 
Zeeman splitting with increasing magnetic field. This 
can be clearly seen in the left panel of Fig. [3] where the 
Zeeman energy, defined as AEz = E(a^ et ) — E(a^ et ) = 
g e x^BB, is plotted for the three dots, as well as for many 
other dots. Positive, negative as well as zero excitonic 
Zeeman splittings were observed for different dots from 
the same sample. We show below that the observed vari- 
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FIG. 2. (color online). Polarization resolved photolumi- 
nescence spectra of three representative quantum dots grown 
without the partially covered island flushing technique. 

ations reflect the range of size and material composition 
fluctuations measured in our X-STM microscopy studies. 
A striking feature in the Zeeman splittings presented in 
the left panel of Fig. [3] is the non-linear dependence of 
AEz on the magnetic field - a magnetic field- dependent 
exciton g-factor g ex = g e -\- g h = g® x + g\ x B. The best fit 
to the Zeeman splittings of QDa, QDb and QDc was ob- 
tained using quadratic function, as depicted by the solid 
lines presented in the left panel of Fig. [3] The solid lines 
for the other QDs interpolate the experimental data. We 
attribute the pronounced non-linear dependence of AEz 
on B to a magnetic field-dependent hole g-factor g^, as 
will be explained in detail below. 

While the excitonic Zeeman splittings of the QDs from 
the non-flushed sample revealed positive, negative and 
even zero B-fleld dependent excitonic g-factors, the dots 
investigated from the flushed sample mainly revealed 
negative, but also B-field dependent g-factors. This is 
clearly demonstrated by the negative and non-linear Zee- 
man splittings presented in the right panel of Fig. [3j The 
weaker variation in the Zeeman splittings observed for 
the dots from the flushed sample compared to the non- 
flushed sample arises from the reduced fluctuation in the 
dot height and the lower average In-concentration caused 
by the PCI growth process. 
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FIG. 3. (color online). Zeeman splitting of the neutral exci- 
ton in several QDs grown without (left) and with the partially 
covered island flushing technique (right). 



A weak non-linear dependence of the electron Zee- 
man splitting on the magnetic field was first observed for 
GaAs-AlGaAs quantum wells and superlattices at high 
magnetic fields P2 Similar quadratic magnetic field de- 
pendencies of the hole Zeeman splitting have been found 
for GaAs-InGaAs superlattices and quantum wells and 
attribute d to a m agnetic field-induced mixing of HH and 
LH states.^ 24 * 25 * In quantum wells subject to magnetic 
field applied parallel to the quantization axis, the hole 
wave function can be factorized in a product of quan- 
tum well states in the vertical direction and Landau lev- 
els in the lateral direction! 25 * 26 * A magnetic field applied 
along the [001] direction induces a coupling of the HH t 
ground state with the LH ^ and the LH \. bands with a 
strength that is proportional to jsy/Bk z and (72 +73) £>, 
respectively.^ The operator k z acts in the growth di- 
rection and 72 and 73 are the Luttinger parameters.^ 
Interestingly, the HH ^ ground state couples only to the 
LH t band via an interaction that varies as (72 — 73) B. 
Since 72 ~ 73 for In^Gai-^As alloys the latter cou- 
pling is normally negligible and the HH | ground state 
consequently has an almost pure HH ^ character, inde- 
pendent of the magnetic field. It has been shown that the 
coupling of the HH t ground state to the light hole bands 
leads to a heavy-hole g-factor that varies with the square 
of the in-plane wave vector, i.e. gu oc /c 2 . 25 For quantum 

well Landau levels, & 2 varies linearly with B leading to 
the experimentally observed quadratic Zeeman splitting. 
On the other hand, in small, strongly confined quantum 
dots & 2 varies with ~ 1/D 2 where D is the dot diameter. 
As a result, the hole g-factor in strongly confined dots 
should, therefore, be unaffected by the magnetic field. 
Our results indicate that the large QDs with dilute In- 
composition produce effects that fall in a regime between 
the expectations for quantum wells (g^ = g® + g\B) and 
quantum dots (gh independent of B). 

To understand the microscopic origin of the observed 



non-linear Zeeman splitting we performed a detailed 
three dimensional electronic structure calculation using 
the eight-band k • p envelope function approximation. In 
order to include the B-field in our calculations we used 
the recently proposed gauge invariant symmetry adapted 
finite element method that accurately accounts for va- 
lence band couplings. ^ Strain fields were included us- 
ing continuum elasticity theory and their impact on the 
electronic structure was fully taken into account via de- 
formation potentials and the linear piezoelectric effect.^ 
The exchange interaction is expected to be of minor im- 
portance for the neutral exciton because of the large ef- 
fective band gap of ~ 1320 meV and the weak mixing 
of conduction- and valence-bands. The direct Coulomb 
interaction was found to have a negligible influence on 
the exciton g-factor and is, therefore, also neglected in 
our simulations.^ 

To obtain quantitative results for the X° g-factor, 
a Luttinger-like eight band k • p- model was employed, 
where remote-band contributions to the effective mass 
Hamiltonian and g-factors are included up to the order 
^,2 127J mo deled our QDs as having a truncated lens 
shape with a diameter varying from D = 25 — 50 nm, 
a height of 4 nm above the wetting layer (WL) and an 
inverse trumpet-like In-compositional profiled The In- 
concentration of the In^Gai-^As alloy was taken to be 
x min _ q 2 a t the base and side of the dot increasing to 
x apex _ 3 _ 0.5 at the dot apexPH These parameters 
are consistent with the results of X-STM measurements 
(Fig. [I]) performed on samples grown under the same 
conditions, from which we also determined the thick- 
ness and In-content of the wetting layer to be 2 nm and 
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FIG. 4. (color online), (left) Calculated Zeeman spin split- 
ting of the neutral exciton in three QDs with the same height 
of 6 nm but different diameter and In-concentration at the 
apex x apex as a function of the magnetic field: (a) 25 nm di- 
ameter and x apex =0.50, (b) 45 nm diameter and x apex =0A5, 
(c) 50 nm diameter and x apex =0A0. (right) Calculated heavy 
hole character of the lowest hole orbital levels having spin-up 
HH t and spin-down HH 4- character, respectively. 

The left panel of Fig. [4] shows the calculated exci- 
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ton Zeeman spin splitting as a function of the magnetic 
field for three model dots having different size and In- 
composition consistent with the range obtained from our 
X-STM measurements. These representative QDs have 
been chosen in order to reproduce the generic behavior of 
QDa, QDb, and QDc, representing the range of behav- 
iors observed in our experiments. The In-concentration 
decreases from the data marked (a) to (c) in Fig. [4] 
from x apex =0.50 to x apex =0A0, whilst the lateral size 
increases from 25 nm (a) to 50 nm (c) to reproduce 
the experimentally observed exciton transition energies 
(1310 — 1365 meV). The curve labeled (a) in the left panel 
of Fig. [4] shows an almost purely linear Zeeman split- 
ting (B-field independent g-factor) whilst (c) exhibits a 
clear quadratic dependence. The model QD (b) shows an 
intermediate behavior between the linear and quadratic 
regimes. The quadratic dependence of the exciton Zee- 
man splitting stems entirely from the HH-like lowest en- 
ergy orbital state in the valence band. For all QDs pre- 
sented in Fig. [4] (a) to (c), the valence band Zeeman split- 
ting varies quadratically with magnetic field, namely, 

AE h z = » B g° h B + » B glB 2 (1) 

where the B-field is applied along the growth direction, 
Hb is the Bohr magneton and, and g\ are the linear 
and quadratic components of the hole g-factor, respec- 
tively. 

These results indicate that the strong quadratic char- 
acter of the hole Zeeman splitting {g\ = 0.016±0.002 T" 1 
for QD A , 0.028 ± 0.002 T" 1 for QD B and 0.026 ± 
0.001 T- 1 for QD C from Fig. || arises from the combi- 
nation of the comparatively large diameter, small height 
and dilute In-content in the present dots. Firstly, the 
dot diameter is larger than the magnetic length over the 
entire range of B-fields of interest. Secondly, the low, al- 
most homogeneous In-concentration induces only a weak 
confinement potential in the core of the QDs. As a re- 



sult, the quantum states that are formed in magnetic 
field resemble 2D Landau levels and, consequently, the 
HH ground states behave in a manner similar to quantum 
wells. In addition, the small dot height of 6 nm (includ- 
ing the wetting layer) introduces a strong field-induced 
coupling of HH ^-LH j~ as in narrow quantum wells 
The quantum- well- like dependence of the Zeeman spin 
splitting on magnetic field is especially pronounced for 
the large, In-dilute QD (c) as shown in Fig. [4] Due to 
the QD shape, In-Ga alloy profile and inhomogeneous 
strain fields, the lowest energy HH orbital has a weak 
LH admixture at zero magnetic field. This is illustrated 
quantitatively in the right panel of Fig. [4] However, as 
the B-field increases, the LH- admixture of the HH ^-like 
ground state increases due to the field-induced LH-HH 
mixing alluded to above. In strong contrast, the HH |- 
like ground state effectively decouples from the LH-bands 
and, thus, its LH character weakens with the magnetic 
field. 

In summary, strongly magnetic field-dependent exci- 
ton g-factors were observed in InGaAs self- assembled 
QDs. The microscopic origin of non-linear Zeeman split- 
ting was accounted for by eight-band k • p simulations 
using realistic parameters (size and In-composition) that 
were directly extracted from X-STM measurements. The 
combined effect of dilute In-composition and relatively 
large dot lateral size was shown to result in strong 
field-induced mixing of the HH-LH orbital states in 
high magnetic fields. This mixing manifests itself as a 
quadratic variation of the hole Zeeman splitting on the 
external magnetic field. Similar effects are negligible for 
the electron and have previously been observed only in 
thin two dimensional systems. 
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